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We investigate theoretically the magnetic levels and optical properties of zigzag- and armchair- 
edged hexagonal graphene quantum dots (GQDs) utilizing the tight-binding method. A new bound 
edge state at zero energy appears for the zigzag GQDs in the absence of a magnetic field. The 
magnetic levels of GQDs exhibit a Hofstadter-butterfiy spectrum and approach the Landau levels 
of two-dimensional graphene as the magnetic field increases. The optical properties are tuned by 
the size, the type of the edge, and the external magnetic field. 
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Graphene is a single atomic layer consisting of a two- 
dimensional honeycomb lattice of carbon atoms. This 
novel system has attracted intense attention because of 
new fundamental physics and promising applications in 
nanoelectronicsiii^ . It exhibits high crystal quality, an ex- 
otic Dirac-type spectrum, and ballistic transport proper- 
ties on a submicron scale. Graphene samples are usually 
fabricated by micromechnical cleavage of graphite and 
have excellent mechanical properties that make it pos- 
sible to sustain huge electric currents. The lateral con- 
finement of Dirac fermions in graphene is still an enig- 
matic and extremely challenging task due to the well- 
known Klein paradox. The Klein paradox makes it im- 
possible to localize the carriers in a confined region uti- 
lizing an electrostatic gate. The confinement of Dirac 
fermions at a nanometer scale is one of the central goals of 
graphene-based electronics and has attracted increasing 
interest ^^i^i^iii^. Recently it was demonstrated experi- 
mentally that graphene can be cut in the desired shape 
and size^'^. Recent progresses in fabricating and charac- 
terizing stable graphene nanostructures provides the op- 
portunity to explore the various remarkable optical^i^iii 
and transport properties^ of these structures. 

In this work, we investigate theoretically the electronic 
structure and optical properties of zigzag- and armchair- 
edged hexagonal graphene quantum dots (GQDs) (see 
Fig. [T]) utilizing the nearest-neighbor tight-binding 
model. The dangling bonds at the edges are passivated 
by hydrogen atoms. The model has been successfully 
used for fullerene molecules, carbon naotubes, and other 
carbon-related maieiiol^i^^^. The Hamiltonian of 
GQDs can be written as = ^e^cjci + Yl ^id^l^j^ 

where Si is the site energy, tij is the transfer energy be- 
tween the nearest-neighbor sites, and c] (q) is the cre- 
ation (annihilation) operator of the tt electron at the site 
i. When considering a magnetic field B applied perpen- 
dicularly to the plane of a GQD, the transfer integral 



tij becomes tij 



te 
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where ^i^- = f J^' • A is 
the Peierls phase. A = (0, Bx^ 0) is the vector poten- 
tial corresponding to the magnetic field B along the z 



axis, which is perpendicular to the graphene plane. In 
our calculation, we take = h/e as the unit of the 
magnetic fiux and <l> = \^BaQ/2 as the magnetic fiux 
through a plaquette, where ao = 2A6A is the lattice con- 
stant of graphite. The difference between the values of 
Si and tij for the atoms at the edge and the center is 
neglected. The relevant parameters used in our calcu- 
lation are e = 0, t = —3.033^^. The eigenvalues and 
eigenstates can be obtained from the secular equation 
det \s-H\= 0, where Ha = 0, H^ij^ = te^^^^^'^' . 

Figs. [H show electronic density distributions of 
the zigzag and armchair-edged graphene quantum dots 
(ZGQD and AGQD, respectively), in the absence of a 
magnetic field. The size of a dot is characterized by N, 
the number of hexagonal units along an edge. Figure 1(a) 
and 1(b) show the probability distributions of the highest 
valence level (HVL) and lowest conduction level (LCL) 
for the ZGQD with smah size {N = 2). The probabil- 
ity distributions of the HVL and LCL correspond to the 
bonding and anti-bonding states that are localized at the 
corner of the hexagonal GQD. In contrast to conventional 
semiconductor quantum dots where the ground state is 
localized at the center of dot, the ground states for the 
conduction and valence bands, i.e., HVL and LCL, lo- 
calize at the middle of each edge in the ZGQD (see Fig. 
[Tfc) and (d)) as the size of the ZGQD increases. This fea- 
ture can be understood as follows: the Dirac fermion in a 
ZGQD behaves like a confined photon in a cavity, and the 
lowest mode is the whispering gallery mode, which also 
localizes at the boundary of the cavity. The difference 
between the bonding and anti-bonding states becomes 
smaller as the size of the ZGQD increases. The differ- 
ence between the edge states of a ZGQD and a graphene 
nanoribbon is that the edge state of the ZGQD localizes 
at the middle of the edge of GQD, in contrast to the ho- 
mogeneously distributed edge state of a zigzag graphene 
nanoribbon^^'^^ or a zigzag triangular GQD^^. This oc- 
curs because the contribution of each carbon atom at the 
edge of a ZGQD to the edge state is different, while it 
is the same for a zigzag nanoribbon or a zigzag triangu- 
lar GQD. The density distributions of the LCL and HVL 
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FIG. 1: (Color online) Electronic density distributions of 
the highest valence level (HVL) and lowest conduction level 
(LCL) in the absence of magnetic field. Panels (a) and (b): 
HVL and LCL for the Nz — 2 zigzag-edged graphene quantum 
dot (ZGQD). Panels (c) and (d) : the same with Nz = 12. 
Panels (e) and (f): the same with armchair-edged graphene 
quantum dot (AGQD). 



in an AGQD (see Figs. [TJe) and (f)) extend more com- 
pletely over the whole GQD region and are very different 
from that in a ZGQD. This difference is indeed caused 
by the different topological geometry of the boundary of 
the graphene nanostructures. 

Figs. [2ja) and (b) show the density of states (DOS) 
of ZGQDs and AGQDs, respectively, with different sizes 
in the absence of a magnetic field. The total number 
of the carbon atoms in ZGQD and AGQD are ^N^ and 
6(3A^^ — ?>Na + 1), respectively. From the figures, we find 
that there is no edge state in a small ZGQD, and the edge 
state appears when the size of ZGQD increases according 
to the states at zero energy. Meanwhile, there is never an 
edge state for the AGQD. To demonstrate how the edge 
state appears, we plot the energy gap, i.e., the energy dif- 
ference between the lowest conduction band level (LCL) 
and the highest valence band level (HVL), as a function 
of the size {N) of the GQD in Fig. [2jc). The energy gap 
decreases as the size of the GQD increases. Interestingly, 
the energy gap of the zigzag (armchair) GQD decays to 
zero quickly (slowly) as the size of the GQD increases. 
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FIG. 2: (Color online) Density of states of ZGQD (a) and 
AGQD (b). We use a Gaussian function f{E) = g-^^-^o^'/^^ 
with a broadening factor F = 0.05 eV to smooth the discon- 
tinuous energy spectra, (c) The gap of ZGQD and AGQD as 
a function of the size, and the function of the fit line \s a/N 
with a 4.9 eV. 



When the size of the AGQD approaches infinity, the gap 
decreases to zero, i.e., we recover the two-dimensional 
graphene case. The calculated energy gap for the AGQD 
falls off as oc 1/L (see the solid line in Fig. [2]^c)), 
where L is the length of each edge of the hexagonal GQD. 
This dependence of the band gap on the size of GQD is 
very different from that of a conventional semiconductor 
QD, which behaves as 

Figs. [3](a) and [D^a) depict the magnetic field depen- 
dence of the energy spectrum of a ZGQD and AGQD ex- 
hibiting a clear Hofstadter butterfly characteristic, which 
is fractal and exhibits self- similar ityi^i^Si^i^^. As the 
magnetic flux increases, the magnetic levels in the GQD, 
i.e., the so-called Fock-Darwin levels, approach the Lan- 
dau levels (see the red lines in Figs. [Sjb) and [H^b)) 
in graphene En = sgn (n) {\^tao/2lB)\/2 |n|, where 
Ib = \lhleB is the cyclotron radius, n is an integer, 
and sgn is the sign function. 

Figs. WLs) [41(c) show in detail how the magnetic 
levels of ZGQD and AGQD approach the zero-th Landau 
level at small magnetic flux. As magnetic flux increases, 
more energy levels approach the zero-th Landau level in 
pairs. The degeneracy of the energy level at zero en- 
ergy will reach its maximum value 2A^ for <l>/$o = 1/2 A^. 
When the magnetic flux increases further, the degeneracy 
of the energy level at zero energy is lifted fast. This fea- 
ture can also be seen in Fig. [3](d), which plots the DOS at 
the Dirac point. This figure indicates that the degener- 
acy, i.e., the number of energy levels at the zero energy, 
approximately decreases inverse linearly with the mag- 
netic flux ^/^o. These figures clearly demonstrate that 
the energy spectrum of the GQD possesses electron-hole 
symmetry when we neglect the second-nearest-neighbor 
interaction. The DOS and the magnetic level fan of the 
AGQD are similar to that of the ZGQD except at small 
magnetic flux. Comparing Fig. [31(c) to Fig. [11(c), the 
magnetic levels in the AGQD are distinct from those in 
the ZGQD at small magnetic flux, because the ZGQD 
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FIG. 3: (Color online) (a) The spectrum of the = 12 
ZGQD in a magnetic field. We use a Gauss function with a 
broadening factor of 0.1 eV to smooth discontinuous energy 
spectra, (b) and (c) the magnetic energy level fan near the 
Dirac point, i.e., the zero energy point. The red lines in (b) 
correspond to the Landau level of two-dimensional graphene. 
(d) the DOS at the Dirac point as a function of the inverse flux 
$/$o, where we use a Gauss function with a small broadening 
factor of 0.01 meV. 



shows the edge state and AGQD does not for the levels 
near the Dirac point in absense of magnetic field. There- 
fore, the magnetic levels exhibit distinct behavior as the 
magnetic flux increases. The DOS of the AGQD (see Fig. 
ID^d)) also shows a step-like feature as the magnetic flux 
at the Dirac point increases. 

Fig. [5] describes the density distributions of the LCL 
and HVL in the ZGQD and AGQD at small magnetic flux 
$/<l>o = 0.01. Interestingly, the density distributions of 
the LCL and HVL penetrate into the center of the GQD 
for the ZGQDs, which is very different from the AGQD 
case where both the electron and hole are dominantly 
localized in the center of the GQD. The density distri- 
butions for the ZGQD and AGQD show Cq^ symmetry. 
This characteristic is caused by the magnetic confinement 
when the magnetic length Ib becomes comparable with 
the size of the GQD. In addition to those differences, the 
LCL and HVL of the zigzag GQD show opposite symme- 
try order with respect to that of the armchair GQD, i.e., 
the LCL (HVL) and HVL(LCL) of the ZGQD (AGQD) 
belong to the Ei{E2) and E2{Ei) representations at zero 
magnetic field (see Fig. [6j). 

The optical properties of GQDs are promising for 
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FIG. 4: (Color online) The same as Fig.3, but for the Na = 9 
AGQD 



potential applications in optic-electronic devices based 
on graphene. Therefore, we calculate the absorp- 
tion spectra of GQD a{f}w) = ^^J^f^^y El^ • Pcv\^ x 

° c,v 

5 {Ec — Ey — hco) ^ where n is the refractive index, c the 
speed of light in vacuum, eq the permittivity of vac- 
uum, mo the free-electron mass, and e is the polariza- 
tion vector of the incident light along the x direction. 
The coupling between the sp2 states and the Pz state 
is neglected since we are only interested in the opti- 
cal properties of the GQD near the Dirac point, i.e., 
at the low energy regime. The momentum matrix^^. 
is {n\p\m) = imo/h^^c*Cr' {r' -r) {pz,r\H\pz,r'). 

r r' 

The momentum operator Px{Py) has E2 symmetry and 
its direct product with all the irreducible representa- 
tions of the Cqv group can be found in Table [H We 
divide the levels of the GQD into two different families: 
Ai,As,Ei e and ^2,^4, £^2 ^ ^2- The symmetry 
requires that only transitions between the valence band 
levels and the conduction band levels belonging to the 
different families Qi and Q2 are allowed. Notice that the 
initial or final states of the transition should belong to 
the El or E2 representations. In Fig. [6Ka) and (d), we 
label the level structure of a A/'^ = 12 and Na = 9 GQD 
near the Dirac point as Ci — Cn for conduction bands 
with ascending order and Vi — Vn for valence bands with 
descending order, respectively. The conduction band lev- 
els Ci and valence band levels Vi belong to the distinct 
families Q.i and 0.2-, respectively. For example, if Ci be- 
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FIG. 5: (Color online) (a) and (b) show the density distri- 
butions of the HVL and LCL for the = 12 ZGQD in the 
presence of the magnetic flex $/$o = 0.01, respectively, (c) 
and (d) the same as (a) and (b), but now for A^^ = 9 AGQD. 



TABLE I: Direct products of the E2 representation for mo- 
mentum operator Px{Py) with all the irreducible representa- 
tions of the Cqv group. The results are presented as direct 
sums of all possible irreducible representations of the Cev 
group. The notations of symmetries are adopted from Ref. 
[23 



Direct product 



Direct sum 



^2(8)^2 
^20^3 
^20^4 

^20^2 



E2 

El 
E2 
El 

^20^40^2 
^10^30^1 



longs to the family l^i, i.e., Ai, As or Ei, Vi must belong 
to the family 1^2, i-e., A2, A4 or E2, or vice versa. For 
zigzag GQDs with even A^^, the conduction band levels, 
from bottom to top, exhibit different symmetries, i.e.. 
As, El, A2 ■ ■ - , the corresponding valence band lev- 
els show El, A4, £^2, Ai For zigzag GQDs with 
odd Nz, the conduction band levels display the opposite 
(same) symmetries £^1, A4, £^2, ^1, • • • to the conduction 
(valence) band levels of zigzag GQDs with even Nz. For 
armchair GQDs, the lowest conduction band level always 
shows the symmetries Ei, A/i, A2, E2^ • • • from bottom 
to top and this order is independent of the size (Na) of 
the armchair GQD. 

For Nz = 12 ZGQD, the lowest optical-absorption 
peak (peak A) corresponds to the transition between 
the lowest conduction band level Ci with E2 symmetry 
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FIG. 6: (Color online) (a) and (d) are the level diagram for 
= 12 ZGQD and Na = 9 AGQD without the magnetic 
field, where different symmetries are represented by differ- 
ent colors and lines: black solid, black dashed, red solid, red 
dashed, green solid, and green dashed lines for the Ei, E2, Ai, 
A2, A3, and A4 irreducible representations of the Cqv sym- 
metry, respectively, (b) and (c) the JDOS and the optical 
absorption spectrum a for A^^ = 12 ZGQD. We used a Gauss 
function with different broadening factors: 0.02 and 0.005 eV 
for the black and red line, (e) and (f) are the same as (b) and 
(c), but for Na = 9 AGQD. 
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FIG. 7: (Color online) The contour plot of the magneto- 
optical spectra of zigzag (a) and armchair (b) GQD, respec- 
tively. 



and the highest valence band level Vi with Ei symme- 
try. The second and third lowest transitions correspond 
to the transition between the level C2 (Cs) with As {Ei) 
symmetry and the level Vs {V2) with E2 {A4) symmetry 
and the level Ci (C4) and V4 (Fi), respectively. But the 
strengths of these three transitions are very small, there- 
fore these transitions are not clearly seen in the contour 
spectrum in Fig. [71 at <l>/<l>o = 0. The strong absorption 
peak (peak D) appears at £^ = 0.26 eV, corresponding 
to the transition between the level Cs with Ei symmetry 
and the level V3 with E2 symmetry. This strong ab- 
sorption arises from the large moment matrix {n\p\m) 
between these states. 
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For Na = 9 AGQD, the lowest peak (peak A) is similar 
to that in the zigzag GQD, corresponding to the transi- 
tion between Ci and Vi. But the second peak (peak B) 
is different from those of the zigzag GQD. This peak cor- 
responds to the transition between the level C2 {C4) with 
A4 (£^2) symmetry and the level V4 {V2) with Ei (A3) 
symmetry. The third strong peak (peak C) indicates the 
transition between the level C4 with E2 symmetry and 
the level V4 with Ei symmetry. Strong absorption takes 
place when the initial (Vi) and final states (Ci) have ei- 
ther El symmetry or E2 symmetry. As the size of the 
GQD increases, the absorption peaks shift to long wave- 
length for both ZGQD and AGQD. The absorption peaks 
of the ZGQD shift to the long wavelength faster than 
those of the AGQD. The relative strength between the 
peak D and A increases as the size of the GQD increases 
for ZGQDs. But for AGQDs, the relative strength be- 
tween the peak C and A is almost independent of the 
size. 

Next, we discuss the effect of a magnetic field on the 
optical spectrum of a GQD. Here, we only focus on the 
small magnetic fiux case (see Fig. [7|). The spectra of two 
distinct GQDs, zigzag and armchair GQD, exhibit quite 
different behavior due to their different level structures 
and the oscillator strengths determined by the boundary, 
especially for the LCL and HVL which localize at the 
edge of ZGQD. The spectra of two distinct GQDs show 
that the strengths of the transitions vary as the mag- 
netic field increases. In particular, the strong absorption 



lines exhibit asymptotic behavior corresponding to 
the transitions between the conduction and valence band 
Landau levels at high magnetic field. We also find anti- 
crossings in the spectra, since the magnetic field induces 
the mixing of the levels belonging to the different fami- 
lies. 

In summary, we investigated theoretically the magnetic 
levels and the optical spectrum in GQDs. In contrast to 
conventional semiconductor QDs, the LCL and HVL ex- 
hibit an edge-state feature, i.e., a non-zero probability of 
being at the edge of the sample, and the density distri- 
bution depends sensitively on the type of boundary of 
GQDs and the magnetic field strength. The magnetic 
levels of GQD display a Hoftstadter butterfiy character- 
istic, and approach the Landau levels of two-dimensional 
graphene as the magnetic field increases. The magneto- 
optical spectrum of a graphene quantum dot in the in- 
teresting energy range (0-3 eV) is promising for carbon- 
based electronics applications. The position and strength 
of the absorption peaks can be tuned by the size of the 
GQD, the type of the edge of the GQD, and the external 
magnetic field. 
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